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Abstract 

A model relating radiative seesaw and minimal dark matter mass scales without beyond the standard 
model (SM) gauge symmetry (Rz^MDM) is constructed. In addition to the SM particles, the Ri/MDM 
contains, a Majorana fermion multiplet Nr and scalar multiplet x that transform respectively as (1, 5, 0) 
and (1,6,-1/2) under the SM gauge group SU{3) C x SU(2) L x U(1) Y . The neutral component N% 
plays the role of dark matter with a mass in the range of 9 to 10 TeV. This scale also sets the lower limit 
for the scale for the heavy degrees of freedom in Nr and x which generate light neutrino masses through 
the radiative seesaw mechanism. The model predicts an iV^-nucleus scattering cross section that would 
be accessible with future dark matter direct detection searches as well as observable effects in present 
and searches for charged lepton flavor violating processes, such as k — > lj"f and u — e conversion. 
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I. INTRODUCTION 



To account for the cosmological observation that our universe is composed of about [H, 0] 
20% dark matter (DM), the standard model (SM) of electroweak and strong interactions has to 
be extended beyond its minimal form. In particular, new particles playing the role of DM have 
to be introduced. One of the most popular candidates is a weakly interacting massive particle 
(WIMP) that leaves a thermal relic density below the relevant freeze-out temperature in the early 
universe. There are many possible ways to introduce such a WIMP candidate. To guarantee that 
the WIMP is stable, most of the models employ a symmetry beyond the SM gauge symmetries, 
such as R-parity in SUSY theories or a Z 2 symmetry in dark matter models with a real or 
complex singlet [if. 

It is interesting to ask whether the SM gauge symmetry alone can already stabilize a new 
particle such that it may be a WIMP candidate. Indeed, it has been shown that such possibility 
can be implemented using the minimal dark matter (MDM) idea jij. The MDM scenario achieves 
stabilization of a WIMP by choosing representations of the SU(2)l x U(l)y gauge group such 
that (a) the choice of hypercharge implies the existence of an electrically neutral component of 
the multiplet, (b) the neutral component has the lowest mass after accounting for radiatively- 
induced mass splittings, and (c) they do not couple to the SM decay products directly at the 
renormalizable level 1 . Therefore, the lightest component field does not decay into SM particles, 
making it a DM candidate. Such a model has an additional bonus that the DM mass is completely 
fixed by the observed DM relic density because the DM thermal relic density is produced by 
known electroweak interactions. 



Apart from the dark matter problem, the observation of neutrino oscillations [1[ necessitates 
extending the minimal SM. The simplest extension introduces gauge-singlet right-handed neu- 
trinos as needed to construct Dirac mass operators. However the vast hierarchy between the 
tiny Yukawa couplings associated with the Dirac mass operators and those associated with the 
charged lepton masses is widely considered theoretically unattractive, and many models, includ- 
ing seesaw models [HH3] and the Zee model jsj, have been developed seeking a more "natural" 
explanation of the neutrino mass scale. Common to these approaches is the existence of new 
particles with large masses above electroweak scale. We refer these large mass scales collectively 
as the neutrino mass scale, which is not determined a priori. 

It has been shown in Ref . that it is possible to link the DM mass scale to the neutrino mass 
scale by generating neutrino masses radiatively, leading to the radiative seesaw model (RSSM). 



This idea has been further studied by various authors [10[. Most of the models proposed impose 
discrete symmetries beyond the SM gauge symmetry to stabilize DM. In this work, we explore 
the possibility of relating DM and neutrino mass scales without imposing additional discrete 
symmetries on the originally renormalizable Lagrangian by combining the MDM and RSSM 
ideas. We refer to this model as the "Ri/MDM" scenario. In particular, we demonstrate that it 
is indeed possible to construct consistent theoretical models to achieve this goal and that they 
may have testable phenomenological predictions. In the following we describe a minimal model 



1 Higher-dimension, non-renormalizablc operators that would generate decays can be suppressed if the associated 
mass scale is sufficiently large. 
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of this type and study some of its implications for charged lepton flavor violation (LFV). We 
show that while this version of the R^MDM is unlikely to be testable with the Large Hadron 
Collider, it can lead to observable effects in future dark matter direct detection experiments and 
searches for LFV. 



II. THE MODEL 



As we explain below, the minimal model of this type contains a new right-handed Majorana 
fermion N R and scalar \ multiplet that transform under the SM gauge group SU (3)c x SU (2) L x 
U(1) Y as 

N R : (1,5,0), x- (1,6,-1/2). (1) 

In addition to the terms already existing in the minimal SM, the most general renormalizable 
interactions containing the Nr and \ are 

L new = N R iYD^N R + (D^xVD^x ~ (L L Y X N R + l -N R MN R + h.c.) - V new 

V new = /4x f x + \«(X^XA + X a Hx (H^H X j x) a + [W#*) 2 + h - c ] , (2) 

where _D M is the gauge covariant derivative and H is the SM Higgs doublet transforming as 
(1,2, 1/2). Here, we include three families of N R but have suppressed the generation indices for 
notational simplicity. The meaning of the a index will be explained later. In the present context, 
we do not consider possible effects of CP-violation, so will take the new couplings to be real. We 
defer a consideration of possible CP-violation to future work. 

If /i^ and the quartic couplings in V new are positive, x wm n °t develop a non-zero vacuum 
expectation value (vev). It is in principle possible that renormalization group running could lead 
to instabilities ((%) ^ 0) at high scales, since x couples to fermions in the third term of L new and 
since fermion loops can cause the quartic scalar couplings to run negative at large field values 



(see Refs. |llHl4| and references therein). However, a suitable choice of the cut-off of the theory 
and tree-level couplings A" and \% x can preclude instability of the (x) = vacuum. The last 
operator in V new , in contrast, receives no such potentially destabilizing contributions at one-loop, 
so the input value of \r x can be chosen at will. As we discuss below, this freedom is critical to 
the viability of the radiative seesaw mechanism. 

The whole Lagrangian is invariant under a Z2 symmetry even after spontaneous symmetry 
breaking ((H) 7^ 0), under which x ~^ ~Xi Nr — > —Nr, and all the SM fields do not change 
sign. Note that this symmetry is not imposed but rather emerges as a consequence of gauge 
invariance given the field content. This accidental Z 2 symmetry makes the lightest particle of 
the component fields in Nr and X stable, thereby providing a WIMP candidate in each multiplet. 
In addition, the active neutrinos do not have masses at tree-level. They arise, rather, at one-loop 
order when integrating out the Nr and x particles. 

The choice of the field content is arrived by the following considerations. The field Nr should 
have zero hypercharge in order to possess a Majorana mass term. One should not allow a 
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LHNr type of Yukawa coupling in order to forbid Nr from mixing with the left-handed lepton 
doublet and thereby allowing the Nr to be stable. This prevents Nr from being a singlet or 
triplet under SU{2) L . The next possibilities are (2,0) or (4,0), where we now drop the SU(3)c 
quantum numbers. However, the component fields in both cases are all fractionally charged 
(±1/2 and ±3/2), thereby precluding the possibility of a DM candidate. In addition, if the 
(2, 0) representation is chosen for Nr, the accompanying Higgs representation x needed for the 
Yukawa coupling term L^Nrx required for radiative neutrino mass generation would be required 
to be (1, -1/2) or (3, -1/2). If (4,0) is chosen for N R instead, \ can be (3, -1/2) and (5, -1/2). 
Again the component fields in x are all fractionally charged. Thus if the lightest particle is in x, 
it will be stable with fractional charge. These models are not phenomenologically acceptable as 
dark matter models. The minimal choice for Nr, therefore, is (5, 0). Its neutral component may 
be the candidate for DM. Once the quantum numbers of the Nr are fixed, the choice for x can 
be made. In order to have the heavy neutrino Yukawa term L^xNr, x must be either (4, —1/2) 
or (6,-1/2). If one chooses (4, —1/2), a WHHx term is then allowed in the Higgs potential. 
This should be forbidden because this term will induce vev for x leading to the breaking of the 
accidental Z 2 discrete symmetry after H develops vev and will also induce x or -^Vr to decay so 
that none of them are stable. 

The fields Nr and x can be written in the tensor forms, Nr = N^^i and x = Xijkim- Here 
the sub- indices take the values 1 and 2, and the fields are totally symmetric under exchange the 
sub-indices. The independent fields (N ++ , N + , N° , N~ , N ) in Nr, and the independent fields 

(x ++ , X + , X°, X~, X~~,X~~~) hi x can be expressed as 

N ini = N ++ , jV 1112 = _LjV+ , N U22 = ^=N°, 

Ni 222 = -^=N~ , N 2222 = N— , 
V4 

1 1 

++ + o 

Xlllll — X > X1X112 — ~7E-X > X11122 — ~~i=X ) 

V5 v 10 

1 _ 1 __ ___ n , 

X11222 — —j=X 5 X12222 — —?eX 1 X22222 — X \ 6 ) 
V 10 v 5 

Writing the Yukawa coupling, the Majorana mass and Higgs potential terms in the tensor 
notation, we have 

L LX N R = L t N jHm x irk >vrn>e j3 'e kk 'e u, e mm ' , N R N R = N^N^e 11 ' e"' e kk ' e 11 ' , 

XX ] = XijklmXljklm > {H ] H X ] X)l = HlxijklmHaXljklm > {H^HX^Xh = H\xijklmHaXl jk lm » 

IT ,.t Ji'a'a 1 tt^.\2 tt tt ,ij ,i' j'kk'll'mm'nn' ( a\ 

H a'Xi'jklm e e ) \ H X) —■HiXjMmn-Hi'Xj'k'l'm'n'CC C e e € , [A) 

where e 12 = 1, e 21 = — 1 and all other elements equal to zero. In the above repeated indices are 
contracted by g l i with g 11 = g 22 = 1 and g 12 = g 21 = 0. 

There are 9 different ways to write (XXX X )<* depending on how the indices are contracted 
by the tensors e v and g l K Using the identity e u e = g lk g^ 1 — g g-* k ', one can show that only 3 of 
them are independent. The independent terms can be chosen to be the following ones 

{XXX X )l = XijklmXijklmXi'fk'l'm'Xi'j'k'l'm 1 > (.XXX V)2 = XijklmX*ijklm'Xi'j'k'l'm'Xi'j'k'l'm 5 
(XXX X )3 = XijklmXijkl'm'Xi'j'k'l'm'Xi'j'k'lm • (5) 
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There are 9 different ways to write (xXX^X^)a depending on how the indices are contracted 
by the tensors e 1 - 7 and g l K Using the identity e^e kl = g lk g^ 1 — g ' g^ k , one can show that only 3 of 
them are independent. The independent terms can be chosen to be the following ones 

(XXX X )l = XijklmXijklmXi'j'k'l'm'Xi'j'k'l'm' > (.XXX X = XijklmXijklm'Xi' 'j'k'l'm'X*' j'k'l'm > 
(XXX X )3 = XijklmX*jkl'm'Xi' j'k'l'm'X*' j'k'lm ' ( 6 ) 

The Majorana mass term for N R expanded in component fields is given by 

N C R MN R = N+ +C MN R ~ - N^ C MN R + N R C MN R - N R C MN+ + N R ~ C MN+ + 

= N—MP R N— + N~MP R N- + N°MP R N° + N + MP R N + + N ++ MP R N ++ (7) 

where Pr denotes the right-handed projection operator and where we have defined the four 
component fields N in terms of Nr by 

N~ = N R ~ + N+ +c , N~ = Nr — N+ c , N° = N° R + N° R C , 
N + =N R C -N R ', N ++ = iV++ + N R ~ C , (8) 

and where M is a matrix in the space of the three N R generations. 

We will work in the basis where M = diag(mjv 1 , m^, ^at 3 ). At the tree level all components 
of Nj (j = 1,2,3) have the same mass m^.. At one-loop level, this degeneracy is lifted with a 
mass splitting between components of different electric charges Q and Q' 0] 

2 

m%-m% = 2 (Q 2 - Q' 2 ) [sin 2 6 w f(m z /m N ) + f(m w /m N ) - f(m z /m N )] , 

f(x) = | |2x 3 In x - 2x + (x 2 - 4) 1/2 (x 2 + 2) ln[(x 2 - 2 - xy/x 2 - 4) /2] } . (9) 

For a given generation, N° has the smallest mass, making it a dark matter candidate. The 
Q — ± partners of Q = component are heavier by about 166 MeV for m N ^> m z making the 
decay N^ 1 — >• A r °7r ± possible. 

After electroweak symmetry breaking in which only the neutral component of H acquires a 
vev, the mass terms for the component fields of x are given by 



m 



lil + \ l Hy v 2 + A 2 ^ 2 



,2 

x— _ ^X ^X 



2,o , 6~ 



2 2,\1 2 , ^ \2 2 , u \ 2 

m fe x o = »x + A ^x t ' + 5 A ^x w + 5 . 

2 6 ~ 

mj mx0 = /i 2 + A^ x w 2 + -A^w 2 - -\h x v 2 , (10) 

where w = (H) is the vev of the Higgs doublet. 

The Xh\ operator mixes x + an d (x~Y > X ++ an d (x Y ■ m the basis of (x )^) an d 
(x ) )i we have the following mass matrices respectively, 



M 



-^\ H y ^ + A^ 2 + fA 2 ^ 2 



I ^x + A V 2 7E~ X "xV 2 

++> ~~ 2 \ 2 ..2 , \1 „2 i 4\2 „,2 ■ V 11 ) 

-^X H% V V X +*H X V +5 A H X V 
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One can adjust the tree level parameters to make the particles above either heavier or lighter 
than the component fields in Nr, and also to split the mass degeneracy of the component fields 
in x- Similar to the case for Nr, loop effects also split the mass degeneracy in the x multiplet. 
For our purposes, the details of the x spectrum are not essential, as we will take x to be heavier 
than Nr (see below). 



III. DARK MATTER 



As mentioned earlier, the choice > implies that x cannot develop a non-zero vev. This 
leads to a residual Z 2 symmetry, under which Nr — > —Nr and x — > — X an d all the other fields do 
not change signs. As a result the lightest component field in x or Nr will be stable and therefore 
can play the role of DM. We emphasize that this residual Z 2 symmetry is a consequence of gauge 
invariance, renormalizability, and concavity of the scalar potential for x and does not result from 
setting any otherwise allowed couplings to zero by hand. 

In general, one could consider either x or Nr dark matter in this model, though for reasons 
we now discuss the Nr case may be more likely to lead to LFV signatures. Since the x field has a 
non-zero hypercharge, its spin-independent direct-detection cross section is governed by tree-level 
Z-exchange. Assuming Q x saturates the relic density, the resulting direct detection cross section 
is too large to be phenomenologically viable 0]. An exception occurs when the mass-splitting 
between the real and imaginar y comp onents of the neutral field, S = Rex° and A = Imx°, 
is greater than about 100 keV 15l4l8l]. In this case, the dark matter particle can only scatter 
inelastically from the nuclei in the detector since the the ZSS and ZAA couplings vanish in 
contrast to the ZSA interaction. For S = \ms —tha] ^ 100 keV, the inelastic scattering process 
mediated by tree-level Z-exchange is kinematically suppressed, yielding a phenomenologically 
viable direct detection cross section associated with the exchange of two gauge-bosons and/or 
the Standard Model Higgs. On the other hand, in order for the S or A to saturate the relic 
density, its mass should be be order lOTeV 19J or larger. These two considerations then imply 
that the coupling |A# X | must be larger than ~ 10~ 4 since from Eq. (ITDl . one has 5 ~ ~ |A# X | — . 

As we will see in Section IV, the product of A# x w 2 and Y 2 sets the scale of neutrino masses, so 
that under these conditions \Y\ < 10 -4 , implying that contributions from this model to lepton 
flavor violating processes will be too small to be observed. 

For a considerably larger value of 5, a lower mass for xo can also be consistent with the dark 
matter relic density and direct detection limits. In this case, the scalar exchange and four-scalar 
interaction contributions to the annihilation process associated with the X^ x terms must be 
taken into account. Cancellations can occur within the allowed parameter space between these 
contributions and those mediated by gauge bosons, leading to the correct DM relic density . 
Detailed studies have been carried out for the inert doublet model. It has been shown that there 
exists two additional mass ranges consistent with the relic density and direct detection limits: 
one with mass below 50 GeV |4|, |20[ and another between mw and 160 GeV 21|. The same 



mechanism can also be applied to the model we study here. However, as with the inelastic dark 
matter scenario, an even larger |A# X | is needed, resulting in even smaller effects in LFV processes. 
Consequently, we defer a study of these lighter-mass x° scenarios to future work, where we will 
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also consider the associated LHC phenomenology. 

Consequently, we will take N R as the possible WIMP DM assuming that the relic density is 
thermally produced in the standard ACDM model. For a consistent model, at least two copies 
of Nr are needed. Only the lightest N R is stable jo), since the heavier N R j can decay through 
the Yuakawa interaction by N R j — > L Li x, followed by x decays into a SM light lepton and the 
N R . If x mass is larger than that of the N R j as we assume here, there is a suppression factor 
due to the off-shellness of x- However, the next-to-lightest N R is still unstable on cosmologically 
relevant timescales. 

To produce the observed DM relic density Vtohih 2 , we require that the thermally averaged 
annihilation rate (ctav) to be 3 x 10 cm 3 s^/flnMh 2 . In our model, a^v is determined by co- 
annihilation 2 of a pair of N into a pair of W gauge bosons and other combinations as shown in 
Fig. [TJ The interactions of Nr with gauge bosons are contained in the kinetic term, Nrj^D^Nr. 
Expanded in component forms, they are given by 

AC NNV = (eA^ + g cos 9 W Z V ){2N++ 1 »N++ + N£-fN% ~ NrI 1 " N r - 2N r ~ ^ N r ~) 
+g(V2N R - + ^W+N R - + y/3N+-fW+N% 

+V3N° R rW+N R + V2N R YW+N R - + h.c), (12) 

which can be written, in terms of the four- component fields N, as 

AC NNV = (eA At + g cos e w Z„){2N ++ -fN ++ + N+-fN + ) 

+g(-V2N ++ YW+N + - V3N + ^W+N + h.c.) , (13) 

where again we suppress the generational indices for simplicity. Note that photon and Z couplings 
to N fields are vector-like. 

Using the above interaction, to the leading order in t> 2 , one obtains jij 

1 2N 1035 4 , 1215 , 
o'av = 7r(c s + c v v ) , with c s = q , and c v = q . (14) 

c s and c p represent the strength of the S-wave and P-wave annihilation. 

By fitting the relic density with the standard thermal relic density calculation, the DM mass 
can be determined. Detailed calculations have been carried out in in Ref.Q, so we will not 
repeat them here. If one requires that Vt DM h 2 = 0.110 ± 0.006 [l| to be saturated by thermal 
annihilation of Nr, the mass of the lightest component must have a mass of 4.4 ±0.1 TeV. A 
more detailed analysis has been carried out in Ref.ji]. Taking into account the Sommerfeld effect 
that enhances the annihilation cross section, the DM mass is raised by about a factor of two 
with |i| rriN = (9.6 ±0.2) TeV. The model has some testable predictions for direct DM detection. 
DM interacts with quarks at one-loop level. The direct DM detection cross section is predicted 
to be ji] 10~ 44 cm 2 which is safely below the current upper limits from CDMSII and XenonlOO 



experiments, but can be tested at future 22[ superCDMS and xenon-lton experiments. 



2 We use the term "co-annihilation" to denote processes in which the N° and either an TV or one of the charged 
components of the multiplct annihilate into SM final states. 
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(a) (b) 




(c) (d) 

FIG. 1: Feynman diagrams for co-annihilations. V and V', / and f are appropriate SM gauge bosons 
and fermions. N% j are component fields in Nr with appropriate combinations of i and j. For example 
the DM component is N{ = N° and its anti-particle Nj = iV 0c . A pair of DM annihilate themselves 
only through the t and u channels in the figure. 
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(a) 

FIG. 2: QfjMh 2 vs. with Sommerfeld effect taken into account. 

The DM mass of order 10 TeV makes it impossible to be directly probed at the LHC. However, 
nature may choose to have several components of dark matter, with N only produce a fraction 
of the total relic DM density 3 . If so, the mass can be smaller. To illustrate, we plot in 
FigfJl the relic density as a function of m^. The relic density contribution from this model scales 
approximately as m 2 N . Therefore, the relic density drops rapidly with decreasing mjv- We see 
that with ttin ~ 1 TeV, the relic density is already only about 1% of the total. If the Ri/MDM 
scenario is going to play a significant role providing the DM relic density, it is unlikely that the 



3 In this case, the model needs to be further extended to accommodate the total DM relic density. 
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Xk > >Xr 
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FIG. 3: One loop Feynman diagram for neutrino mass generation. 
LHC will be able to probe the heavy degrees of freedom our particular model realization. 

IV. NEUTRINO MASSES 

Light neutrinos do not have masses at tree-level because the vev of x is zero. However, at one- 
loop level non-zero neutrino masses can be generated [s| through the diagram shown in Fig. |3j 
The relevant operators are the Nr Yukawa interaction and the \h x {Hx) 2 + h.c. operator. None 
of the other quartic interactions involving H and \ that appear in V new can lead to a lepton- 
number violating neutrino mass operator at one-loop since they contain a x and X* pair. The 
corresponding vertices in the diagram can be obtained from the following terms 

LlxNr u\^=Nffx" + \J\ N R 3 X~ + \J\n%X° ~ ^ N RjX + + A^X ++ ) 



{H X f -+ - A -^X + X~ + fxY) , 



(15) 



where the indices i and j are generation indices. 

There are three pairs which can contribute to light neutrino masses in the loop. They are 
(N ++ , N ), (N + , N~), and (N 0c , N°) with appropriate pairs of x component fields in the loop. 
The light neutrino mass matrix is given by 

, . 1 14 YljYlH ~ ry . 

{M v )iV = ~rp, o ~r *h x v I{m k ,m m ,m Ni ) , 

16n 2 5 m Ni 

1 i m X^ m Xra^ Nl ) = — '-j-l-XO-z ^- 2 , 16 

mi — mi \ mi T — mi mi r — mi J 

Xk Xm \ JVj Xk JVj Xm ' 
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where rn Xkm and m^. are the masses of the scalars, and fermions in the loop. 

Neglecting the mass splitting of the component fields in x (with a common mass m x ), 
I(mk,m m ,m N ) becomes 

/(2) = a-g^in, (17) 

where z = m 2 N /m x . In the limit = m x , 1(1) = 1/2. 

Diagonalizing M u , one obtains the Pontecorvo— Maki— Nakagawa— Sakata (PMNS) mixing 
matrix Vpmns f° r the charged current interaction involving leptons, M u = Vpmns^v^rmns- Here 
M v is the diagonalized light neutrino mass matrix. The signs of the eigen-masses can be nor- 
malized to be positive by appropriate chiral rotations of the heavy neutrino fields. A general 
solution for Y can be written as 

/ 7, 2 14 \ " 1/2 - 

Y = ( ^i~*H X ) V VM ^Ml /2 OM]l 2 I^I\ (18) 

where M n and / are diagonal matrices with = diag(m^ l ,m^ 2 , ...) and / = 

diag(I(zi), I{z2), ■■■)■ O satisfies 00 T = I and can be complex in general. 

From the above one finds that in order to have at least two massive light neutrinos, one 
requires at least two copies of Nr fields. With two Nr, one of the light neutrino masses is zero. 
If it turns out that all the three light neutrinos are massive, at least three Nr fields must be 
introduced. In the presence of multiple, non-degenerate Nr generations, the lightest one plays 
the dominant role in setting the mass scale for neutrinos. If they are almost degenerate, all will 
contribute significantly. In all these cases, the DM mass plays the role in setting the lower scale 
of the heavy particles in the RSSM. 

With three degenerate Tn^ 12 ^ = m^, one can easily find solutions of Y which give the correct 
neutrino mixing pattern and also their masses in the allowed ranges. As an illustration, we 




take = 1 and \ Hx = 10" , m N = m x = 9.6 TeV, Vpmns to be the tri-bimaximal form [23 
and obtain sample Yukawa couplings 1nh,ih with the light neutrinos in a normal mass hierarchy 
(0,0.009, 0.05) eV and an inverted hierarchy (0.05,0.051, 0.01) eV. We have the corresponding 
Yukawa couplings given by 

0.032 \ / —0.11 0.078 

Y NH = \ 0.032 0.095 , Y m = 0.055 0.078 0.046 (19) 
0.032 -0.095 J \ 0.055 0.078 

For different choices of O, the resulting Y are also different, reflecting the fact that the 
number of parameters is larger than the number of constraints. The key point, however, is that 
for a generic choice of O one can have large Yukawa couplings that may in turn have some 
implications for other leptonic flavor violating phenomena. Importantly, the Yukawa coupling Y 

7—1/2 ~ 

scales as X H ^ . By adjusting the size of \h X i one can change the overall magnitude of the Yukawa 
couplings. In the above a small number for A# x has been used. If one sets this coupling to be 
zero, the model has a global U(l) lepton number symmetry. In that sense a small number for 
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Xh x is technically natural. Indeed, explicit inspection of the one- loop corrections shows that this 
operator does not mix with the others in V new under renormalization. Moreover, as emphasized 
earlier, the a non- negative Xh x wm remain so after RG running. So the scenario with small A h x 
is consistent with vacuum stability considerations. 

As discussed above, for \Xh x \ having a significantly larger magnitude (of order ~ 10~ 4 ), one 
may also be able to consider inelastic x dark matter, since in this case the associated mass- 
splitting between the real and imaginary parts of x° will be sufficiently large to suppress the 
inelastic direct detection cross section. However, from Eq. ( Fl8l) and the overall scale of neutrino 
masses, we find that the magnitude of the Yukawa couplings must be generically < 10~ 4 . As we 
will see in Section |V] the corresponding effects on lepton flavor violating observables would then 
be too small to be observable in the next generation of LFV experiments. Since our emphasis 
here falls on the possible signatures of the R^MDM scenario and LFV, we will not elaborate 
further on the large inelastic \ dark matter possibility. 



V. LFV: PRESENT CONSTRAINTS AND FUTURE PROBES 



With the dark matter mass determined 10 TeV, this model has zero chance to be tested by 
directly producing the new heavy degrees of freedom at the LHC. As discussed earlier, a N R mass 
lighter than one TeV that could in principle lead to LHC discovery would imply undersaturation 
of the relic density by a factor of 100 - too small to play a significant role for dark matter. As 
an alternative probe, we consider charged lepton flavor violating processes assuming that this 
model saturates the DM relic density. Indeed, the light neutrino masses arise at one-loop order 
and have a magnitude governed by M N and M x and by products of the Yukawa couplings Y and 
\h x , whereas loop-induced LFV amplitudes are independent of \h x , Consequently, it is possible 
that that Yukawa couplings can be sufficiently large to generate observable LFV signatures while 
remaining consistent with the scale of light neutrino masses. In what follows, we show that the 
model can lead to a large \i — > branching ratio, comparable in magnitude to the recent limit 
reported by the MEG collaboration 24| and its expected future sensitivity 251 as well as a large 



26 , COMET 27 , and PRISM 28 



e conversion rate that may be accessible with the Mu2E 
experiments. The LFV decays r — > /i(e)7 can also provide additional information. 

Constraint from fi — > e-f 

The decay amplitudes for ii — > ij^f can in general be parameterized as 



M 



6 £j % (J 



^{ArPr^A^U 



(20) 



Neglecting small mass splitting of component fields within a given Nr multiplet, we obtain 



.4 



X- k 



Y,Y* k Y jk F(z k ) 



F(z) = P } 



2z 3 + 3z 2 - 6z + 1 - 6z 2 In z 
6(1 -zY 



P 



N ' 



6^ 2 + 3^ + 2 + 6^1n2 
6(1 -zY 



(21) 
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where Zf. = m 2 Nk /m x . Summarizing all one loop contributions, we have P x = —5 and P/v = 2. 
Al can be obtained by replacing by rrij in the above. 

Note that the fj, — > e^y amplitude does not depend on the parameter X Hx . As emphasized 
earlier, this feature allows the possibility of having small A# x , but large Y to satisfy constraints 
from neutrino masses and mixing, and to have a large \i — > e^y branching ratio. Moreover, since 
the DM relic density is governed by the gauge rather than Yukawa interactions, one still requires 
knowledge of the to determine Zj — » Ijj, even with the scale of heavy degrees of freedom 
fixed by the DM relic density. It is clear that the known constraints on light neutrino masses 
and mixing matrix provide important information, yet they do not complete fix the products 
Y* k Yjk- Consequently, one may consider the LFV searches as independent probes of the Yukawa 
structure of the model, assuming sufficient sensitivity. To illustrate and to simplify the analysis, 
we will take the heavy degrees of freedom to be degenerate, i.e., = m^.. Note that the mass 
splittings between generations so as to allow for a one-species DM scenario can be sufficiently 
small that we may neglect them here. Similarly, the splittings within a generation induced by 
weak radiative corrections are negligible for this purpose. We will then work with the partial 
branching ratio defined as BR(Zj — > Ijj) = Br(Zj — > Zj7)/Br(Zj — > l^vv). 

Although there presently exists no evidence for any decays k — > lj~f, impressive bounds have 
been obtained in many cases. For many years, the most stringent on BR(/i — > erA is 1.2 x 10 _11 [l|. 
Very recently, the MEG collaboration has obtained a better upper limit with|24| BR(/i — > ej) < 
2.4 x 10 -11 at 90% C.L. One way to ascertain whether our model can have testable consequences 
is to see if with known constraints, the model can produce partial branching ratios close to the 
current experimental bounds. To this end, we first vary m vx from eV to 0.08 eV for the normal 
hierarchy, and 0.05 eV to 0.09 eV for the inverted hierarchy, with the light neutrino mass square 



differences, Am 2 21 = (7.59±0.21) x 10~ 5 eV 2 [3l|, and |Am§ 2 | = (2.43±0.13) x 10" 3 eV 2 [32J fixed 



from experiments. These choices also satisfy the cosmological constraint Y2i m ^i < 0-28 eV|33l|. 
We then determine the allowed parameter space assuming that this model saturates the Br(yU — > 
ej) experimental upper bound 2.4 x 10~ 12 . We then obtain the constraint between |A# X | and 
m x with Eq. ffl8|) applied as well. The results are shown in Fig. HJ 

We comment that in order to have a large branching ratio for /j, — > e^y close to the upper 
bound in both normal and inverted hierarchies with m Vl larger than 0.08 eV, the parameter \h x 
is smaller than 10~ n . In this regime, the elements of the Yukawa matrix Y are of order a few, 
which may be close to the non-perturbative region wherein the results may become less accurate 
compared with smaller m Ul . But for m Ul less than 0.08 eV, the perturbative predictions should 
be reasonably reliable. 

In Figj5l we also show the constraint on the relevant product of Yukawa couplings using the 
current experimental bounds on /i — > &y and \i — e conversion. Present constraints from fi — > 
as well as searches for \l — > e conversion in nuclei (see below) are given in the left panel, while 
the sensitivities of prospective experiments are shown in the right panel. We see that for m x 
not too much larger than m^, that is z = (mN/m x ) 2 close to 1 the combination of the Yukawa 
coupling J2 k Y e kY* k is presently constrained to be less than 1 (left panel). The upper curve in 
the right panel gives the expected sensitivity of the MEG experiment. We observe that even 
with the assumption that the Yukawa couplings are less than order one, the model may have 
observable effects in ji — )■ ej. 



12 



(a) Normal Hierarchy 



(b) Inverted Hierarchy 



- m v , =0eV 
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»«,. =0.08eV 





0.2 0.4 



FIG. 4: log 10 \ ^Hx\ vs - z = m N/ m x ^ or ( a ) norma l hierarchy, and (b) inverted hierarchy for different 
choices of m Ul . The curves are obtained using the recent MEG upper limit Br(/i — > ej) = 2.4 x 10' 
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FIG. 5: Comparison of present (left panel) and prospective (right panel) constraints from fi—e conversion 
on various nuclei and fj, — > &y. 



Constraints from \x — e conversion 

In the case of \l — e conversion, new operators contribute to the LFV amplitude; in addition to 
the dipole interaction (liZUj) one also must include the charge radius operator. In general there is 
also a Z-penguin contribution. However, due to the vector-like coupling of Z to N, this operator 
vanishes for the same reason that [icy charge operator vanishes, a feature that we have verified 
explicitly. This situation differs, for example, from that in the SM with four generations, where 
Z-penguin dominates the contribution 29|. Moreover, the Z dipole and charge radius-induced 



ji — e conversion amplitude is suppressed by a factor of m^/m|. Consequently, all Z-exchange 
contributions can be safely neglected. 

The photon charge radius-generated quark level [/, — e conversion amplitude given by 

L = E Y ekY; k ej,PMP x G x (z k ) + P N G N (z k )] QqQl^Q , (22) 
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where Q q is the electric charge of quark q in unit e, sw = sin Ow and cw = cos 9w, and 



G x {z) 



G N (z) 



2-9z + 18z 2 - llz 3 + 6z 3 In z 

36(1 -zY ' 
-16 + Ahz - 36z 2 + 7z 3 + 6(3z - 2) In* 
36(1 -z) 4 



(23) 



Several groups have performed searches for \x — e conversion Q. In Tabled we list results for 
the conversion-to-capture ratio 



ni r(/i- + A(N, Z)^e-+ A(N, Z)) 

^ e Ti pt r( Ji -+A(N t Z)^v lt + A(N+l,Z-l) 



(24) 



where A denotes the atomic number 



To obtain /x — > e conversion rates on different nuclei from the FCNC interaction in the above, 
we start with the effective four-fermion efiqq operators that may contribute. Following the 
notation of Ref. 30 we have 



^cff 



4GV 

'7f 



[m^ u (A R P R + A L P L )fiF^ u + h.c. 



G 

[e(9LS( q )PR + 9RS( q )PL)H qq + e(g LP{q) P R + g RP{q) P L )fi qj 6 q + h.c] 

G 

-j= [Kqlv^Pl + 9RV{q)l^PR)^qi»q + e(9LA(q)YPL + 9 RA(q)Y p r) V Tl m7s2 + h.c] 



Gp 
72 



^e(g LT(q) a flu P R + g RT ( q )^ 1 'Pj)\tqp >g + h.c. 



(25) 



Comparing with Eq. [22j we have, at the one-loop level, g(L,R)(s,p,A,T)(q) = 0, g R v(q) = 0, and 

V2 A R 



A R 



9LV{q) 



8 GfttIh 



m e , 

A L = —A R 



m, 



'u 



2tt 2 



J2^Q q Y ek Y; k z k [P x G x (z k ) + P N G N (z k )} . 



(26) 



(27) 



The contribution from At can be neglected compared with that from A R . 

Theoretical efforts have been made by several group to calculated these matrix elements. We 
will use the results in Ref. 30| for a consistent calculation. The resulting expression for quantity 



P>^e measuring the leptonic FCNC effect in fi — > e conversion is proportional to 



A R D + g®V® + MvV>\ 2 



9lv* 



(28) 



where D(A), V^(A) and V^ n \A) are overlap integrals as a function of atomic number 30J and 

9lv = 2 9lv(u) + gLV(d) , g ( Lv = 9lv(u) + 2g LV{d) . (29) 
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A 


upper limit 


D(A) 


t r( rA / A \ 

V (P) (A) 


t r/nl'l / A \ 

y m, (A) 


<^ e ( A ) 


27 a 1 
13 A1 




0.0362 


0.0161 


0.0173 


0.0026 


32 q 


7.0 x 10" 11 


0.0524 


0.0236 


0.0236 


0.0028 


|Ti 


4.3 x 1(T 12 


0.0864 


0.0396 


0.0468 


0.0041 


79 97 AU 


7.0 x 10~ 13 


0.189 


0.0974 


0.146 


0.0039 




4.6 x 10~ n 


0.161 


0.0834 


0.128 


0.0027 



TABLE I: The relevant parameters for /i — e conversion processes, which are evaluated by using method 
I in Ref.^ 



301 ] . The upper limits are their 90% c.l. values. 



Combining these, we obtain 
Bi 



ArD(A) A r D(A) 



(30) 



R°(A) = L-Jt—\D(A)\ 2 . (31) 



B(fJ. ery) 

where 



The relevant parameters are listed in Table [J, which are evaluated using method I in Ref . |30| • 

The constraints on Ylk ^fc^ujfc f rom At — e conversion on various nuclei and \i — \ cy are shown 
in Fig|S] on the left panel assuming that are degenerate. We see that among the current 
available experimental limits, fi — > cy gives the strongest constraints. The current bounds 
on fi — e conversion do not provide significant constrains on the parameter since the Yukawa 
couplings can still be much larger than of order 1. On the other hand, the sensitivities of 
prospective future fi — e conversion searches Mu2E[26| /COMET [2~7| and PRISM[28j exceed that 
of the MEG experiment. 

Effects on r — > /i(e)7 

We now study the allowed ranges for r — > fi(e)y. The current experimental bounds are: 
BR(r — > /i7) < 2.5 x 10~ 7 and BR(r — »■ erf) < 1.8 x 10~ 7 [if. For illustration we will consider the 
case for degenerate Ni with = 9.6 TeV. In this case once the value of Br(/i — > e*y) is fixed, 
the entries in the Yukawa matrix Y are also fixed in terms of m Ul . Therefore, predictions for 
Br(r — > fiy) and Br(r — > ej) can be made as a function of m Ul . We show the results in Fig. [6] and 
Fig. [71 From Fig. [61 we see that Br(r — > /17) increases, and decreases when m Ul increases for the 
normal and inverted hierarchy cases, respectively. The predicted partial branc hing ratio can be 



about 10 which may be probed by future experiment at the super B factory [341] . For r — > ey 
the partial branching ratio is predicted to be much smaller than the current experimental limit 
making it difficult to test the model using this process. 

In the examples given above, we have used the tri-bimaximal form for the mixing matrix 
Vpmns- We have also checked cases with 6*13 set to be non-zero. We find a similar pattern, but 
slightly smaller partial branching ratios. Because in our model, we have approximately, BR(Z,; — > 
Iji) = BR(/i -> ey)\[{VpMNsYM u (VpMNs)%\V\[(VpMNsVM u (VpMNsy}i2\ 2 , with anon-zero l3 
the denominator of the right-hand side becomes larger. This leads to smaller branching ratios 
for r —> /i(e)7 compared with the cases discussed earlier. 
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FIG. 6: Br(r — > ^7) as a function of m Vl in (a) normal hierarchy, and (b) inverted hierarchy. 
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FIG. 7: Br(r — > ej) as a function of m Ul in (a) normal hierarchy, and (b) inverted hierarchy. 



Flavor diagonal lepton magnetic moments 

Before ending this section, we would like to make some comments on magnetic dipole moments 
of fi and neutrinos. The anomalous magnetic dipole moment of muon a M can be generated at 
one loop level in our model. The current value of difference between experiment and SM is [l[ 
Aa^ = (255 ± 63 ± 49) x 10 -11 . The new contribution for a/ =ejAtjT in this model is given by 
Aag = (mj/167r 2 m 2 N )zF(z) J2k l^fc| 2 - Again, for a given m ui with fixed fi — > er) branching ratio, 
Aa,£ can be obtained. The result is much smaller than the current constraint given previously. For 
example, Aa^ = — 6.8x 1CT 11 when m Ul = 0.08 eV with normal hierarchy, and Aa M = — 1.5x 10~ n 
when m Ul =0.05 eV with inverted hierarchy. 

Neutrinos can also have magnetic dipole moments. For Majorana neutrinos only transition 
magnetic moments ^ may be non-vanishing. In our model, the one loop level calculation, 
obtained by attaching an external photon to the charged lines in FigJ3] in all possible ways, 
results in a vanishing /ijj. It can only be generated at higher orders. We have estimated the two 
loop contribution and found that is of order 10 -20 fiB- This is much smaller than current 



experimental upper bound is [i v < 0.32 x 10 Hb at CL = 90% [35 
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VI. SUMMARY 



In this work, we have constructed a scenario - the Rz/MDM model - relating the "minimal" 
dark matter and radiative seesaw neutrino mass scales. An important feature is that this model 
does not impose any beyond-the-SM gauge symmetry. It contains, in additional to the SM 
particles, a Majorana fermion multiplet Nr and a scalar multiplet x which transform respectively 
as (1,5,0) and (1,6, -1/2) under the SM gauge group SU(3) C * SU{2) L x U{1) Y - The lightest 
new component of the Nr is electrically neutral and does not decay into SM particles, making it 
a natural candidate for DM. The DM can (co-)annihilate into SM particles through electroweak 
interactions. To produce the correct DM thermal relic density, the DM mass is determined to be 
in the range 9 to 10 TeV. This scale also sets the lower limit for the scalar x mass scale which, 
in combination with Nr, generates light neutrino masses through radiative seesaw mechanism. 
The large DM mass of order 10 TeV makes it impossible to directly detect the new particles at 
the LHC. This model, however, predicts that there is a sizable direct DM detection cross section 
of order 10~ 44 cm 2 which can be tested in future experiments, such as superCDMS, and Xenon- 
lton. This scenario can also accommodate present neutrino data while leading to a large \i — > e'y 
branching ratio that can be probed by the MEG experiment. Present constraints obtained from 
fx — e conversion searches are weaker than those from fx — > e'y. However, future fx — e conversion 
experiments on Al and Ti can provide much more stringent constraints. This model also can 
imply a sizeable r — > fi'y decay rate that may be observable at super B factories. 
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